
sigma = diag(1, 2, 2)
mu1 = c(-1, 0)
mu2 = c(1, 0)
mu3 = c(0, 1.7)
N = 30